Tunnelling density of states at Coulomb blockade peaks 



O 
O 



43 



a 



o 



X 



Nicholas Sedlmayr, Igor V. Yurkevich and Igor V. Lerner 
School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, United Kingdom 

(Dated: 25 July 2006) 

We calculate the tunnelling density of states (TDoS) for a quantum dot in the Coulomb blockade 
regime, using a functional integral representation with allowing correctly for the charge quantisation. 
We show that in addition to the well-known gap in the TDoS in the Coulomb-blockade valleys, there 
is a suppression of the TDoS at the peaks. We show that such a suppression is necessary in order to 
get the correct result for the peak of the differential conductance through an almost close quantum 
dot. 



*0 ! I. INTRODUCTION 



The Coulomb blockade in quantum dots is one of the most thoroughly investigated and best understood phenomena 
in modern condensed matter physics (for reviews seei2). Starting from the earliest papers on the subject )?) 4 )^? , the 
attention was focused on the conductance through Coulomb-blockaded dots and on the corresponding statistics of 
heights, widths, and positions of the Coulomb-blockade peaks in conductance. A complimentary approach was to 
describe the tunnelling density of states (TDoS) in the c\o&£&. While also leading to the description of conductance 
and its distribution it was of considerable interest by itself. The reason is that a singularity in the TDoS in the 
zero-dimensional dot potentially bridges the zero-bias anomaly in the TDoS of interacting electrons in the metallic 
regime 10 with the Coulomb gap in the insulating regimeAi. 

Qualitatively, behaviour of the TDoS found inLSai has a clear origin. The Coulomb blockade makes the energetic cost 
of putting an extra electron in the dot of the order of the charging energy, E c — e 2 /C (C is an effective capacitance). 
' Consequently, it leads to a low-energy gap in the TDoS stretching up to E c at low temperatures, T <C E c . By tuning 
the gate voltage applied to the dot one periodically reaches degeneracy points where the energy of having N and 
N + 1 electrons in the dot is the same. At these points, i.e. at the peaks in the Coulomb blockade regime, tunnelling 
O [ is no longer suppressed and one expects the gap in the TDoS to close. In this Letter we show that, although the 
TDoS at the peaks remains finite, its zero-energy value is of half of its value at e > E c . We then demonstrate that 
y—i \ such a suppression is necessary to restore the correct results&^ii for the peak conductance. 

J> . We will use technique similar to that introduced by Kamenev and Gefen- in their original calculation of the TDoS 
0^ ' in the quantum dot, representing the electron Green's functions as functional integrals, albeit we use the Keldysh 
representation in the form developed wil rather than the Matsubara one as irA The main difference in our approach 
is that in the saddle-point approximation, we take into account all the electron winding numbers. This is necessary 
to get the correct expression for the TDoS at the Coulomb blockade peaks and also is the only consistent way (in 
the discussed context) to get the electron number quantization both at the peaks and in the valleys. Naturally, this 
could be obtained either in Keldysh or in Matsubara formalism. The necessity of including all the winding numbers 
has already been stressed ioi£ for the problem of a granular medium comprised of an array of quantum dots in the 
Coulomb blockade regime. However, the effective averaging over chemical potentials of the dots inevitable for the 
granular system, gives in this case only the value of the TDoS in the valleys. 

We consider the problem in the ergodic (zero-dimensional) limit. It is straightforward in our formalism to include 
a correction from spatially inhomogeneous modes; however, this correction is small as l/<? (where g ~ 1/t^ w 8 3> 1 
is the dimensionless conductance of the dot and Td w is the electron's mean dwelling time in the dot) and will not 
differ from that found previously^. Thus we start with the standard 'universal' Hamiltonian for a zero-dimensional 
• • ■ system2i2ii4, keeping there only the charging term: 



H = H + ^(N-N g ) 2 . (I) 



Ho is the Hamiltonian of free electrons in a random potential V, N is the electron number operator, eN g is the 
neutralising background charge (governed by the gate voltage for the standard quantum dot). It is convenient to 
represent the appropriate Green's function in terms of the functional integral in the Keldysh technique 1 ^, 

iG(r, t; r', t') = Z' 1 [ V^Vifj i/>(t, t)tp{r' , t') e iS M , (2) 







FIG. 1: The 'interaction' Keldysh contoui^. 



where Z = f Vi^Vi/j e iS ^ and the action S[ip] = S [ip] + S c [ip] is given by 
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SoM = J dtj dr^(r,t) \id t - £j ^(r,t), 
SM=~^J dtN 2 (t), N(t)= J dr^(r,t)^(r,t) 



(3) 



The time integrals are taken along the standard 'interaction' contou»i^i&, FigQ] and Ji = \i + E c N g on the thermo- 
dynamic (imaginary) part of the contour and /i = E c N g on its dynamical (real) part. This corresponds to a Green's 
function defined as an average with the grand canonical Gibbs density matrix; its time evolution is described in the 
Heisenberg representation with a Hamiltonian, Eq. |T|I. that does not contain \x. In the zero-dimensional regime 
considered here it is convenient to expand Grassmann fields ^(r, t) and ip(r, t) in terms of free electron eigenfunctions, 
introducing fields c(t) and c(t) which depend on time only: 

^( r ,*) = ^2^n{r) c n (t), |Vn( r ) =£n^n(r), = e« - M- 

n 

The Green's function becomes 



In the zero-dimensional system under consideration positions r and r' in Eq. are indistinguishable so that observable 
quantities can be found from 



G(t, = /" d d r G(r, i; r, = £ G n (t, t') . 



(5) 



Decoupling the charging term 5* c in Eq. Q with the help of the standard Hubbard-Stratonovich transformation leads 
to (i) replacing S c by the bosonic action 



s[4>] = ~J dt<p{t) 



(6) 



and (ii) substituting idt — <f> for idt in the action Sq. To calculate the functional integral over the fermionic fields c n (t) 
and Cn(t), and thus G(t,t'), we notice that for any function tp(t) 



Z^ 1 I Vc n Vc n c n (t)c n (t')e*-(K dt5 ^ 9t+ ^ c - = sgn(MV J '' dr ^ (r) 



The time ordering here is along the contour in Fig. ^ so that the integral in the second line above is also taken along 
the contour and sgn(t, t') equals 1 (or —1) when t precedes (or goes after) t' on the contour. The bosonic field is not 
included in the normalization factor Zq in the above expressions. 

Note the analogy with the calculations in the Matsubara technique in Ref& the bosonic field can be gauged out 
by a shift in the fermionic field but for the initial conditions in the imaginary time integral. These do not allow one 
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to get rid in this way of the zero-frequency field component. In our calculations <f>o is the precise analogue of this 
component. Hence 

G n (t, f ) = -* s y) / vjeiSM ~ n{M e ti , z = | H(<w . (7) 

Here we have defined the grand canonical partition function, S(0o), and the grand canonical partition function with 
the n th level excluded, E n ((j)o), with energy levels in both shifted by the charging effects: 

S(0 O ) = JJ (1 + e -«»+*>) , S„(0 O ) = S(0 O ) (1 + e-«-+*')- 1 . (8) 

71 1 

It is convenient to expand E(0o) an d £ n (0o) in Eq. (HJ in terms of the canonical partition functions: 



S„(0 O ) = f] Z w (e n ) e^+^ , Z N (e n ) = <f g e" Wfl JJ I 1 + e-*»+«) 
The meaning of Z/v(£„) can be seen from its formal definition: 



(9) 



^v(e„) _ Tr ^ {°»<i e ~ 



Tin (e-^o) 



1-Fjv(e„), (10) 



where F^(e n ) is the distribution function in the system of iV noninteracting electrons, as the charging energy in the 
Hamiltonian is just a constant when the number of electrons N is fixed. Note that in order to get a consistent 
description of the charge quantization in the quantum dot, it is crucial to have all winding numbers m in the expansion 
©• 

On substituting this expansion and a similar one for S n into Eq. J7J), and performing the Gaussian integration over 
the fields <j), we find all the Keldysh components of the Green's function. Thus, the Fourier transform of G > (t,t') 
(where t is on the lower, and t' is on the upper part of the Keldysh contour, Fig^l is given by: 

9 ' °° 

G > (e) = -^J2/2 e- 0EN Z N (e-n N )5(e-e n -n N ), (11) 



Z 

n N=0 



where 



Z = 52 e-fi E » Zir , E N = ^ (N - N g ) 2 - p N , Cl N = E c (JV + | - N g ) 



2 

N=Q 

Now we can effectively average over disorder simply by substituting the mean TDoS of noninteracting electrons, i/q, 
for S(e — e n — fijv) , with the assumption that the TDoS is smooth in any realisation of disorder which is valid 
when the mean (or typical) level spacing 8 is much smaller than T. Next we use Eq. I|1(J|I to obtain 

G>(e) = -^p £ e ~ PEN t 1 ~ F ^ « ~^p- E e ^ EN [1 - f( £ > ( 12 ) 

N N 

where we have used that Fn(s — £In) for N ^> 1 is approximately the same as the (grand canonical) Fermi-Dirac 
distribution function /(e — f2jv) with the 'chemical potential' /l*o of order NS (which is small as compared to f2jv 
and thus discarded). Note that this obvious relation between canonical and grand canonical expressions can be 
easily confirmed directly by using the definition l|l()(l and calculating Zn and ^jv(£ n ), Eq. J^J, in the saddle point 
approximation. 

The TDoS can be found from the standard formula 

v{e) = ^[G R (e) G A (s)} = ^[G>(e) - G<(e)} . (13) 
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(c) 

FIG. 2: The dependence of the TDoS (in the units of vo) on the energy (measured in E c ): (a) in the valley, (b) through an 
intermediate region, and (c) at the peak. 



Although G < and G > are exactly related (in the equilibrium case) by G^(e) 



W G„ (e) , this is not convenient 



to use for approximations in the strong Coulomb blockade regime (E c (3 3> 1) as both functions in the product are 
very sharp. Instead, we note that there exists an exact relation, Zn = Z^{e n ) + e~P En Z^^i{e n ) that follows from 
the definition (|5J), which allows us to express G < after straightforward transformations as 



G<(e) 



E 

N 



-f3E* 



/(e-O^_0. 



From these expressions for the Green's functions we find the TDoS to be: 

' e ~" EN [/ ( £ - E ^ N -h~N g ))+l-f(e- E C (N + I N g ))] . 



N 



Keeping only the leading terms in the sum over N, we get 



v{e) U{e-n N ) + e~^ n «-^U(e-n N+1 ) 



(14) 



(15) 



(16) 



where we have denned U(e — SItv) = /(s — ^iv-i) + 1 — /(e — Qjy). All other terms are exponentially suppressed so 
that we just have the terms for N closest to N g + i. Away from the degeneracy point, one of the terms in eq. (|16l) is 
also exponentially suppressed so that this expression is contributed by one term only. This corresponds to the valley 
in the Coulomb blockade and gives the TDoS with a gap, as illustrated in fig. |2Ja). In approaching the degeneracy 
point, fig. I3b), the gap is smeared by the contribution from the other term in eq. (|16f) : finally, at the degeneracy point 
corresponding to the peak of the Coulomb blockade, the TDoS remains finite at all energies but shows a 'half-gap' at 
e < E c , fig. E[c). 

Now we show that the half-gap in the TDoS is necessary to restore a correct expression for the current through 
an almost closed quantum dot. We consider such a dot connected via weak tunnelling contacts to leads (reservoirs) 
with fixed chemical potentials. Introducing in the usual wajii& the retarded, advanced and Keldysh components of the 
Green's function, one can rewrite the standard expressionii for the current through the a th lead (with r Q = 2irv a \^ a | 2 
being the tunnelling rate, and j a tunnelling coefficients) as follows: 



r + co j 

I a = eT a J Tr [G K (e) - [l - 2f a (e)] [G R (e) - G A (e)] } 



(17) 



Restricting considerations to the case of two leads (a = 1,2) and using the current conservation, X\ + = 0, 
exclude the Keldysh compponent G K from eq. I|17|) to find 



we 



T — 



eTiF 

r 



f + OO J 

j ^-[h^)-fM[G R {e)-G A {e)] 



(18) 



In the linear response regime, when the difference between the chemical potentials in the leads ^2 — ^1 — &V , we 
obtain the following expression for the differential conductance, 



di 2 rir 2 



As 



df(e - /x) 



de 



^0 



(19) 
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Here v(e) is the TDoS calculated for the dot in a thermodynamical contact with the leads at the same chemical 
potential /i. Substituting the expressions for the TDoS obtained above and keeping only the two leading order terms 
N and N + 1 in the TDoS, i.e. using eq. Hlfi|) . we arrive at 



+00 



ds 



df(e - aQ 
de 



v{e) n N /2T 



vq sinh J7jv/2T ' 



provided that |fijv| = \E a (N + |) — N g \ <C E a . Substituting this into eq. (|TT?)l reproduces the standard resuh-^ for 
the differential conductance: 

di 2 r x r 2 n N /2T 

av= e «W2f' (20) 

We note that for reproducing the correct coefficient at the peak in the differential conductance it is necessary to have 
the half-gap structure in the TDoS described above. 

In conclusion, we have calculated the Green's functions describing an isolated quantum dot in the Coulomb blockade 
regime. This allows us to find a complete description of the TDoS both in the Coulomb valleys and peaks. At the 
peaks we find a new feature in the TDoS: the zero-energy suppression of the TDoS by the factor of 2 in comparison 
to its value at high energies. By considering the quantum dot weakly coupled to two leads we can find the linear 
response conductance for elastic tunnelling. We have also demonstrated that the correct description of the TDoS 
at the peak is necessary to obtain the standard result for the differential conductance. Finally, we believe that the 
techniques developed here could be extended for considerations of the spin blockade effects in the quantum do1ji£. 



Acknowledgments 

This work was supported by the EPSRC grant GR/R95432. 



1 ALHASSID Y., Rev. Mod. Phys. 72 (2000) 895. 

2 Aleiner I. L., Brouwer P. W. and Glazman L. I., Phys. Rep. 358 (2002) 309. 

3 Kulik I. O. and Shekhter R. L, Zh. Eksp. Teor. Fiz. 68 (1975) 623. 

4 Averin D. V. and Likharev K. K., J. Low Temp. Phys. 62 (1986) 345. 

5 Mullen K., Gefen Y. and Ben- Jacob E., Physica B 152 (1988) 172. 

6 Glazman L. I. and Matveev K. A. JETP Lett. 51 (1990) 484. 

7 Nazarov Y. V., Zh. Eksp. Teor. Fiz. 96 (1989) 975. 

8 Kamenev A. and Gefen Y., Phys. Rev. B 54 (1996) 5428. 

9 Levitov L. S. and Shytov A. V. JETP Lett. 66 (1997) 214. 

Altshuler B. L. and Aronov A. G., Sol. State Commun. 30 (1979) 115. 

1 Shklovskii B. I. and Efros A. L., J. Phys. -Condens. Matter 8. 

2 Kamenev A. and Andreev A., Phys. Rev. B 60 (1999) 2218. 

3 Efetov K. B. and Tschersich A., Europhys. Lett. 59 (2002) 114; Phys. Rev. B 67 (2003) 174205. 

4 Kurland I. L., Aleiner I. L. and Altshuler B. L., Phys. Rev. B 62 (2000) 14886. 

5 Keldysh L. V., Zh. Eksp. Teor. Fiz. 47 (1964) 1515. 

6 Rammer J. and Smith H., Rev. Mod. Phys. 58 (1986) 323. 

7 Jauho A. -P., Wingreen N. S. and Meir Y., Phys. Rev. B 50 (1994) 5528. 

8 Kiselev M. N. and Gefen Y., Phys. Rev. Lett. 96 (2006) 066805. 



